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Abstract 

We study double Higgs production at the future Linear Collider in the framework of the Two 
Higgs Doublet Models through the following channels: e + e~ — > &i&jZ, <3?j = h°, H°, A , H^. 
All these processes are sensitive to triple Higgs couplings. Hence observations of them provide 
information on the triple Higgs couplings that help reconstructing the scalar potential. We also 
discuss the double Higgs-strahlung e + e~ — > h°h°Z in the decoupling limit where h° mimics the 
SM Higgs boson. The processes e + e~ — ► h°h°Z and e + e~ — > h°H°Z are also discussed in the 
fermiophobic limit where distinctive signatures such as 4j + X , 27 + X and 67 + X are expected 
in the Type-I Two Higgs Doublet Model. 

PACS numbers: 12.60.-t, 13.66.Fg, 13.85.Lg 
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I. INTRODUCTION 

In order to establish the Higgs mechanism for the electroweak symmetry breaking 
we need to measure the Higgs couplings to fermions and to gauge boson as well as the 
self-interaction of Higgs bosons. Such measurements, if precise enough, can be helpful in 
discriminating between models through their sensitivity to quantum corrections, particularly 
in specific scenarios such as the decoupling limit where Higgs couplings mimic the SM ones. 

If electroweak symmetry breaking is achieved by the Higgs mechanism, it is possible to 
discover at least one light Higgs boson at the Large Hadron Collider (LHC). Such a Higgs 
boson can be a Standard Model (SM) one or one of those predicted by various extentions 
of the SM, such as the Minimal Supersymmetric Standard Model (MSSM) or Two Higgs 
Doublet Model (2HDM). 

In case of a discovery of the Higgs boson at the LHC, it may be possible to measure its 
couplings to gauge bosons and fermions at a certain precision [2| . It has been demonstrated 
in Ref. [3|] that physics at the LHC and the e + e~ International Linear Collider (ILC) will 
be complementary to each other in many respects. In many cases, the ILC can significantly 
improve the LHC measurements. 

In recent years there has been growing interest in the study of extended Higgs sectors 
with more than one Higgs doublet 0, 0, The simplest extension is the 2HDM; and such 
a structure is required for the MSSM. Models with two (or more) Higgs doublets predict the 
existence of charged Higgs bosons. Therefore, the discovery of them would be the conclusive 
evidence of an extended Higgs sector. 

For the linear collider, there have been several studies dedicated to triple Higgs couplings 
both in the SM and beyond (for a review, see Ref [3]). In Refs. [8| and jj the double 
Higgs-strahlung processes and WW fusion both in the SM and MSSM have been addressed. 
In both cases, the sizes of double Higgs-strahlung processes and WW fusion cross section 
are rather small. It has been shown in Ref. 10J that for the process of e + e~ — > ZHH with 
H — > bb in the SM, the irreducible background from electroweak and QCD processes can be 
suppressed down to manageable levels by using kinematics cuts. 

At the LHC, the double Higgs production has been studied in Ref. 11| with the conclusion 
that triple Higgs couplings in the SM and MSSM can be measured provided the background 
can be rejected sufficiently well. 
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In 2HDM's, triple and quartic Higgs couplings have been studied for a linear collider 



in Refs. 
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In Ref. 



121 ] . the triple and quartic couplings have been studied in 
2HDM without CP violation. The analytic expressions of triple and quartic couplings are 
given, but numerical evaluations for the cross section of double Higgs-strahlung are given 
in the framework of MSSM. Ref. [13J| studies the triple Higgs couplings in the framework 
of 2HDM with CP violation in the Higgs sector. Besides, the numerical analysis for triple 
Higgs couplings has been presented in MSSM with explicit CP violation. 



Ref. [14( concentrates exclusively on the triple Higgs production of the 2HDM at the 
linear collider, and the cross sections are found to be large. In this paper, we study the 
double Higgs-strahlung production e + e~ — > Q^jZ, = h°, H°, A , in the general 
2HDM, as measurements of these processes can shed some light on the 2HDM triple Higgs 
couplings. In addition, we also take into account the perturbativity as well as vacuum 
stability constraints on various parameters in the Higgs potential. We will show that after 
imposing those constraints, there are still large enough cross sections, a few hundred femto- 
barns (fb) in some cases, to allow a determination of the 2HDM triple Higgs couplings. We 
will also study some of these processes in the decoupling limit and in the fermiophobic limit 
of the so-called type-I 2HDM. 

The paper is organized as follows. In the next section, we give a short review of 2HDM 
as well as a rough numerical estimate of the triple Higgs couplings. In Section IHIt we 
evaluate the double Higgs-strahlung production e + e~ — > &i&jZ in the 2HDM paying special 
attention to e + e~ — *> h°h°Z in the decoupling limit, where h° mimics the SM Higgs boson, 
as well as in the fermiophobic limit. Our findings are summarized in Section HVl 

II. SHORT REVIEW OF 2HDM AND TRIPLE HIGGS COUPLINGS 
A. Short review of 2HDM 



In this section we define the scalar potential to be studied in this article. It has been shown 

n n, 

in Refs. |4|, [6| that the most general 2HDM scalar potential that is both SU(2)l <g> U(l)y 
and CP invariant is given by: 

$ 2 ) = m?$I$i + m^l$ 2 - (m? 2 $!$2 + h.c) + -Ai($l$i) 2 + -A 2 ($£$ 2 ) 2 
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+ A 3 ($ 1$ 2 ) + A 4 ($ t 1 $ 2 )($ T > 2 ) + iA 5 [($i$ 2 ) 2 + h.c] , (1) 

where $1 and $ 2 have weak hypercharge Y — 1 and vacuum expectation values (VEV's) t>i 
and w 2 , respectively, and Aj and m\2 are real- valued parameters. Note that this potential 
violates the discrete symmetry $j — > — $j softly by the dimension-two term m 2 2 ($i~$ 2 ), and 
has the same general structure as the scalar potential in the MSSM. 

After the electroweak symmetry breaking, the W ± and Z gauge bosons acquire their 
masses. Explicitly, three of the eight degrees of freedom in the two Higgs doublets correspond 
to the three Goldstone bosons (G ± , G°) and the remaining five become physical Higgs 
bosons: h°, H° (CP-even), A (CP-odd), and H ± with masses m h o, m H o, itiao, and m#±, 
respectively. 

The potential in Eq. ([1]) has ten parameters (including v\ and v 2 ). The parameters mi 
and m 2 are fixed by the minimization conditions. The combination v 2 = v \ + v\ is fixed as 
usual by the electroweak breaking scale through v 2 = (2\/2(jf) _1 . We are thus left with 
seven independent parameters; namely (\i)i=i t ...,5, rriu, and tan/5 = v-ijv\. Equivalently, we 
can take instead 

m h o , m H o , m^o , m H ± , tan/3 , a and nii 2 . (2) 

as the seven independent parameters. The angle (3 diagonalizes both the CP-odd and charged 
scalar squared-mass matrices and a diagonalizes the CP-even squared-mass matrix. One can 
easily calculate the physical scalar masses and mixing angles from Eq. in terms of Aj, 
77ii2 and Vj. and invert them to obtain Aj in terms of physical scalar masses, tan/3, a, and 



m 12 
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16]. 



It is straightforward to derive the triple Higgs couplings from the above scalar potential in 
Eq. (Cp). In the next section, we list the trilinear scalar self-couplings relevant for our study. 
Other relevant couplings involving Higgs boson interactions with gauge bosons and fermions 



can be found in Refs. [4j, [17|]. We note that Ref. [18| lists the complete Higgs trilinear and 
quartic interactions for two types of 6-parameter potentials, referred to as 'Potential A' and 
'Potential B'. Potential A is equivalent to our potential if m 12 — > 0, and in this limit the 



Feynman rules in the next section are in agreement with those in Ref. 181 ] . 



There exist three classes of 2HDM's. The main difference among them is in the ways they 
couple the Higgs fields to matter fields. Assuming natural flavor conservation [jjj], the two 
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most popular models are type-I and type-II, denoted by 2HDM-I and 2HDM-II, respectively. 
In 2HDM-I, the quarks and leptons couple only to one of the two Higgs doublet exactly the 
same as in the SM. In 2HDM-II, one of the two Higgs fields couples only to down-type 
quarks (and charged leptons) and the other one only couples to up-type quarks (and neutral 
leptons) in order to avoid the problem of large flavor-changing neutral currents (FCNC's). 
This is also the pattern found in the MSSM. The third class of models is type-Ill, denoted 
by 2HDM-III. However, they are generally regarded as problematic as FCNC's appear at 
the tree level [2^]. We note that 2HDM-I can lead to a fermiophobic Higgs boson h° in the 



limit where cos a = 
this model is h° — > 77 or h° 
not have such a feature. 



211 ] . and the dominant decay mode for the lightest Higgs boson in 
WW, depending on its mass. In comparison, 2HDM-II does 



In our analysis we also take into account the following constraints when the independent 
parameters are varied. 



3iggs scalars 22| should not 
23| : |Ap| < 10~ 3 . Such an 



The extra contributions to the Ap parameter from the 
exceed the current limit from precision measurements 
extra contribution to 5p vanishes in the limit m H ± = m^o. To ensure that Ap be 
within the allowed range, we demand only a small splitting between m H ± and m^o. 



From the requirement of perturbativity for the top and bottom Yukawa couplings 24] 
tan/3 is constrained to lie in the range 0.3 < tan/3 < 100. 



We note in passing that it has been shown in Ref. [25|] that the latest B — > Xg'y 
branching ratio among others puts a lower bound on the charged Higgs mass, vtlu± > 
250 GeV, in 2HDM-II. However, the conclusion does not apply to 2HDM-I. Since our 
analysis is generally valid for both 2HDM-I and 2HDM-II, we still consider values of 
mji± below 250 GeV for that matter. 

To constrain the scalar-sector parameters we will use both perturbativity constraints 
on A, 



26 



271 ] as well as vacuum stability conditions 



28 



291 ] . We require that all 



quartic couplings of the scalar potential in Eq. (JTJ) remain perturbative: Aj < 8n for 
all i. These perturbative constraints are slightly less constraining than the full set of 
unitarity constraints 26[ [27I established using the high energy approximation and the 
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equivalence theorem. For vacuum stability conditions, we use 28l. l29j]: 



Ai > , A 2 > , 
^/X^+ A 3 + min(0,A 4 - |A 5 |) > (3) 



From the experimental point of view, the combined null searches from all four CERN 
LEP Collaborations give the lower limit m#± > 78 GeV at 95% confidence level (CL), 



a limit which applies to all models satisfying BR(if ± — » rv T ) + BR(if ± — > cs)=l 30]. 
Two LEP Collaborations (OPAL and DELPHI) have performed a search for a charged 
Higgs decaying to AW* (assuming m^o > 2mb) and derived limits on the charged 



Higgs mass 30| comparable to those obtained from the search for H ± — > cs, rv. 



For the neutral Higgs bosons, the OPAL and D 



on the masses of h° and A in the 2HDM 



31 



LPHI Collaborations have put a limit 



321 ] . OPAL concludes that the regions 



1 < m h o < 55 GeV and 3 < m^o < 63 GeV are excluded at 95% CL, independent of 



a and tan/3 3JJ. DELPHI Collaboration studies the Higgs to Higgs decay h° — > A A 
in e + e~ —>■ h°Z and h°A° production and a large portion of the m^o-m^o plane is 
excluded, depending on the suppression factor that enters the cross section formulas 



321 ] . In what follows, we will assume that all Higgs masses are greater than 100 GeV 



except in the fermiophobic limit where m^o can be as light as 60 GeV. 

In Fig. [U we show the allowed region in the m 12 -tan/3 plane, taking into account the 
theoretical constraints mentioned above as well as Ap. The surviving parameter region is 
quite restricted on the tan f3-rrii2 plane. The perturbativity and vacuum stability constraints 
together dramatically reduce the allowed parameter space of the model. In particular, the 
perturbativity constraint excludes large values of tan/3. In fact, with sin a = 0.6, m h o = 120 
GeV and for m#o = 240 GeV to 480 GeV, the upper bound on tan /3 decreases from 10 to 
5.4 at m 12 = 0. Nevertheless, for specific values of m 12 7^ and m#, one can see that values 
of tan f3 as large as 70 can still survive the perturbativity and vacuum stability constraints. 

However, in the fermiophobic limit of 2HDM-I (right panel of Fig. [1]) where a = tt/2, 
m/jo = 60 GeV and 200 GeV < m#o < 500 GeV, tan f3 can be as large as about 30 for m 12 < 
15 GeV. Note that both plots in Fig. [T]are symmetric under m 12 — > — m 12 . For subsequent 
analyses, we will take tan (3 = 10 as a typical value. 
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FIG. 1: Allowed region on the m^-tan/? plane, taking into account perturbativity and vacuum 
stability constraints. The common parameters are m^± = 250 GeV and m^o = 150 GeV. The 
left panel uses (sin a, myfi) = (0.6, 120 GeV) for differents values of mjjo, and the right panel is for 
200 GeV < m H o < 500 GeV and (sina,m h o) = (1,60 GeV) 




120 140 160 180 200 220 240 260 10 20 30 40 50 60 

m H o (GeV) m I2 (GeV) 



FIG. 2: The 2HDM tree level self couplings \f$£j% normalized to the SM self coupling \fjf h : (left) 
as a function of m H o with m\2 = GeV, (right) and as a function of m\i with nrijjo = 240GeV. 
The others parameters are m^o = 120 GeV, mjj± = 140 GeV, m^o = 150 GeV, tan/3 = 10 and 
sin a = 0.8 
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B. Triple Higgs couplings in 2HDM 



In this section, we will study the behavior of triple Higgs couplings as a function of the 
2HDM parameters m h o, m H o, m A o, m#±, tan/3, a and m\ 2 . These couplings are given by: 
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(4) 
(5) 

(6) 
(7) 
(8) 
(9) 
(10) 

(11) 
(12) 
(13) 

(14) 



where \fplp h o is the SM triple Higgs coupling. Here we use the short-hand notations: sw = 
singly with 8w being the weak mixing angle, and and denote respectively sin0 and 
cos0 for <j) being various linear combinations of a and (3. 

As we can see from Eqs. (!5l)- (fl4l) . all triple Higgs couplings have some quadratic de- 
pendence on the physical masses m$ and soft breaking term m\ 2 . These couplings also 
depend strongly on tan f3 and a. In the present study, we will show that varying the scalar 
parameters within the allowed range can still make the triple Higgs couplings larger than 
the SM triple Higgs coupling Afo^ 0/l o = — 3em^ / (2mwSw) by several orders of magnitude. 
These 2HDM triple Higgs couplings can also be large compared to the MSSM triple Higgs 
couplings. The reason for this is that in the MSSM, supersymmetry imposes restrictions on 
the quartic couplings Aj by relating them to the gauge couplings 6, §}. 
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FIG. 3: The 2HDM tree level self couplings A^^f normalized to the SM self coupling \^{ as a 
function of tan/3. The plots use m h o = 120 GeV, m H o = 370 GeV, m A o = 150 GeV, m H ± = 250 
GeV, m\2 = 46 GeV, sin a = 0.6 (left panel) and sin a = 0.7 (right panel). 

Taking into account all the previous theoretical and experimental constraints listed in the 
above subsection, we give here the sizes of the triple Higgs couplings involved in the double 
Higgs-strahlung production. 

In Fig. [21 we illustrate the sizes of the triple Higgs couplings as a function of CP-even 
Higgs mass m#o (left plot) and as a function of the discrete symmetry breaking parameter 
m\2 (right plot). One can see from these plots that the sizes of triple Higgs couplings 
normalized to the SM triple Higgs coupling lie between —10 and 12 for all allowed values of 
m#o and niyi- 

If the splitting between H ± and A is small, it follows from Eqs. and (fT5T) [Eqs. (TTU|) 
and ([T4l ] the couplings \h°A°A° and X h o H + H - [Xh°a°a and X H o H + H -} are almost equal. The 
couplings Xh0H+H- and Xh°h+h- are thus not shown. It is clear from Eqs. ([5]) and ffTUj) that 
A/iO/jO/jO and A^o^o^o have no m#o dependence. That is why \hPh°hP and XhOA°A° exhibit no 
variation over m#o. For the other couplings, the variation over m#o can change the size of 
these couplings by about one order of magnitude. 

The dependence of the triple Higgs couplings on myi is also important (see the right 
plot of Fig. [2]). However, due to the perturbativity and vacuum stability conditions, m 12 is 
constrained to be less than 55 GeV. It is interesting to note that the couplings Xh0h°h° and 
X H o h o h o contributing to e + e~ —>■ Zh°h° enter with opposite signs and may have a destructive 
interference in the amplitude. The same observation holds for X^aoa , Xhoa°a°, Xhon+H- 
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FIG. 4: Double charged and neutral Higgs-strahlung production in the 2HDM at e + e linear 
colliders with S = h° or if . 

and Xh°h+h-- 

In Fig. [31 we show the variation of triple Higgs couplings as a function of tan j3. We fix 
m 12 = 46 GeV, m h o = 120 GeV, m H o = 370 GeV, m A o = 150 GeV, m H ± = 250 GeV, 
and sin a = 0.6 for left panel and sin a = 0.7 for right panel. Almost all the triple Higgs 
couplings are enhanced in the large tan/9 limit. The triple Higgs couplings X H o H o H o, X h o h o h o, 
^H°h°h° an d Xh°H+H- can be more than 50 times larger than the SM triple coupling. This 
tan/3 effect has been first noticed in Refs. [3^]. It is also worth pointing out that some 
of the 2HDM couplings Xhihjhk can flip sign as one varies myi and/or tan/3. 

Note that only in Fig.[3]we ignore perturbativity and vacuum stability constraints. In fact, 
in the left pannel (right panel), perturbativity constraints are violated for 6 < tan/3 < 49 
(7.5 < tan/3 < 48) while vacuum stability constraints are violated for 43 < tan/3 < 49 
(36 < tan/3 < 50). 
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FIG. 5: cr(e + e~ — > H + H~ Z) in units of fb as a function of center-of-mass energy (left) for several 
values of m H ± and as a function m H ± (right) for ^fs = 800 GeV. Regarding the other parameters, 
the left plot uses m h o = 120 GeV, m H o = 370 GeV, m\2 = 45 GeV, m^o = 150 GeV, tan [5 = 53 
and cos(/3 — a) = 0.5. From top to bottom in the right plot, (tan/3, mi2)=(55,44 GeV), (40,52 
GeV), and (30, 60 GeV). 

III. PROBING TRIPLE COUPLINGS FROM DOUBLE HIGGS-STRAHLUNG 
PROCESSES 



In this section we will cover the following processes: e + e — > H + H Z, e + e — > h°h°Z, 
e+e" -> h°H°Z, e+e" -> , and e+e" -> The two processes e + e~ -> 

and e + e _ — > H°A°Z are not sensitive to any triple Higgs couplings; they proceed only 
through gauge couplings and will not be calculated here. 



Measurements of those processes will give some information on the invo 



couplings. The complete calculation is done with the packages FeynArts 
and Loop Tools 



ved triple Higgs 



34|, FormCalc |35|, 



In these 2 — > 3 processes, a width for the internal Higgs exchange is 



needed to stabilize the phase space integration. Such a width is introduced and computed 



at the tree level. The numerical eva 
done with the help of CUBA library 



uations of the integration over 2^3 phase space is 



37). 



Before discussing our numerical results, it is worth pointing out that the following results 
are valid for both 2HDM-I and 2HDM-II since they involve only the Higgs coupling and the 
gauge coupling to Higgs or to e + e~ . 
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Feynman diagrams for the first process e + e~ — > H + H~ Z are depicted in the first set of 
Fig. HI Only diagram (g^) is sensitive to the triple Higgs couplings \h°H+H- and \r°h+h-- 
The other diagrams proceed through gauge interactions only. Similarly, Feynman diagrams 
for the processes e + e~ — > $>i<&jZ with $3 = h°,H°,A° are depicted in the second set of 
Fig. |H As we can see, only diagram (6 2 ) has a triple Higgs couplings dependence. The other 
diagrams depend only on gauge couplings. 

To illustrate the size of e + e~ — > H + H~Z cross section, we show in Fig. [5] cr(e + e _ — > 
H + H~Z) as a function of center-of-mass energy for mjj± = 100,150,200 GeV. The cross 
section exceeds 100 fb for light charged Higgs mass of 100 GeV. For moderate charged Higgs 
mass of 200 GeV, the cross section is still of the order of a few fb. 

In the right panel of Fig. we illustrate the sensitivity of the cross section to m H ± for 
a fixed center-of-mass energy of 800 GeV. By varying m#± from 100 to 350 GeV, the cross 
section is suppressed by several orders of magnitude. However, for a charged Higgs mass in 
the range 100 to 300 GeV, it is still possible to have a total cross section larger than 0.1 fb. 
This can lead to 100 raw events for the planned luminosity of 1000 fb" 1 . 

The collider signature for the e + e~ — > H + H~ Z process depends on how H ± decays. 
Below the top-bottom threshold and before the charged Higgs to neutral Higgs decays 
H ± —> SW ± (S = h°,A°) are open, the signature would be Zt + t~ + $ T or Zcscs. 
Once the charged Higgs to neutral Higgs decays H ± — > W ± S become open, one would get 
the ZW + W~ SS final states with S decaying either to bb or t + t~. Above the top-bottom 
threshold, the signature would be Ztibb if Br{H ± — > tb) dominates or ZW + W~ SS if 
J Br(i7 ± -> W ± S) dominates. 

Let us now turn to double Higgs-strahlung processes e + e~ — > Qi§jZ, with $jj being 
neutral Higgs bosons h°, H°, and A . As noted in the previous section, the &iH + H~ and 
§iA°A° couplings (<J>j = h°, H°) are almost identical for not too large mass splitting between 
charged and CP-odd Higgs. As we will see, the total cross section of e + e~ —>■ ZA°A° has 
similar behavior and size to the total cross section of e + e~ — > H + H~Z for m^o and m#± of 
about the same magnitude. 

In Fig. [61 we show the total cross sections for the neutral modes Zh°h°, Zh°H°, ZH°H° 
and ZA°A° as a function of the center-of-mass energy yfs for different values of m^o. As 
explained in the figure caption, we select specific sin a values to optimize the total cross 
sections. The behavior of the cross section for those neutral modes is similar to that of 
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FIG. 6: Total tree level cross section of e + e~ — > Z<&i$i in 2HDM as a function of y/s with 
m h o = 120 GeV, m H ± = 140 GeV, m 12 = GeV, tan/3 = 10 for different values of m^o. The left 
plot uses nifjo = 200 GeV and sin a = 0.6. The right plot uses m H o = 240 GeV and sin a = 0.8. 

e + e~ — > ZH + H~. The reasons are that for both processes we have similar topologies 
in the contributing Feynman diagrams and that the triple Higgs couplings have a similar 
magnitude. For the Zh°h° and ZA°A° modes, the total cross section is maximized for 
y/s ~ 350 — 400 GeV. In the case of Zh°h° with m h o = 120 GeV, the maximum of the 
total cross section is reached for m^o = 150 GeV, m H o = 240 GeV, and sin a = 0.8 (right 
panel), and is slightly larger than 100 fb. This enhancement for e + e~ — > Zh°h° is due to the 
threshold effect (H° -> h°h° is open). In the ZA°A° mode, the cross section is maximal for 
m^o = 150 GeV, m#o = 200 GeV, and sin a = 0.6 and is about 300 fb. The cross section 
of the heavy mode ZH°H° is very small and is not shown; whereas the cross section of the 
mixed mode ZH°h° is mild, fa 1 fb, in both cases. 

In order to show the sensitivity of the cross sections to high tan/3, we plot in Fig. [7|the 
cases of small tan/3 (left) and large tan/3 (right). For small tan/3, e + e" — > ZA°A° dominates 
with a few fb cross section, e + e~ — > Zh°h° and e + e~ — > Zh°H° barely reaches 0.1 fb. This 
is due to the fact that the triple Higgs couplings involved in these processes are such that 

(\2HDM \2HDM \2HDM \ 2HDM \ / n O. e n A iq 1 ! v \ SM 
I'VhO/iO, AH°h°h°i A h°A°A°i A H°A°A°) ~ (U.O, 0, U.4, Id) X A h o h o h o ■ 

They are not enhanced by much except for A^o^o^o which is 13Afofo fe o and explains the 
dominance of e + e~ — > ZA°A° in this case. However, in large tan/3 cases, e + e~ — > Zh°h° 
has the largest cross section (fa 60 fb), followed by e + e~ — > ZA°A° (fa 20 fb). In this large 
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FIG. 7: Total tree level cross sections of e + e — > Z&i&i in 2HDM as a function of y/s with 
m h o = 120 GeV, m H o = 370 GeV, m H ± = 250 GeV, m A o = 150 GeV, m 12 = 46 GeV, sin a = 0.6 
and tan/3 = 5 (left) or tan/3 = 50 (right). 

tan p case, the involved triple Higgs couplings are given by 

(\2HDM \2HDM \2HDM \ 2HDM \ ^ f ao on or\\ v \SM 
\ A h,0h°h°i A H0h,0h°i A h°A°A°i A H°A°A°) ~ {^6, —6V, OZ, 6\J ) X A h o h o h o ■ 

Obviously, they are much more enhanced than the low tan f3 case. 

Fig. [H] shows the behaviors of e + e~ — > Zh°h° and e + e~ — > ZA°A° (left plot) and the 
branching ratios of the CP-even H° (right plot) as a function of m H o and for m h o, m H ±, 
m^o, mi2 = 120, 250, 150, GeV, respectively, tan/5 = 5, sin a = 0.9 and y/s = 500 GeV. 
To compute the cross sections, we have summed the continuum part a(e + e~ — > ZSS) and 
resonant part generated in the chain e + e~ — > ZH° — > ZSS via the resonant H° Higgs- 
strahlung, with S = h° or A . It is obvious that the cross sections of e + e~ — > Zh°h° 
and e + e~ — > ZA°A° reach their maxima near the threshold regions of H° — > h°h° and 
H° AM . 

In the right panel of Fig. [HI we show the decay branching ratios of H°. The bb mode 
dominates below the WW threshold, while the WW mode takes over for m H o > 2Myy. 
Above the h°h° threshold, the decay H° — > h°h° is open and dominates over the others at 
the order of 70%. For m#o > 300 GeV, H° — > A A is open and has also a substantial 
branching ratio of (9(20%), reducing the branching ratio of H° — > h°h° to the level of less 
than 20%. Note that once H° — > H + H~ is open, it quickly reaches the level of 15%. 

As shown in the left plot of Fig. [HJ there are three kinks occurring when m#o = 
2m h o,2m,AO,2mH± and corresponding to the opening of the H° — > h°h°, H° — > A A and 
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FIG. 8: Left plot: total cross sections for Zh°h° and ZA°A° continuum and resonant production: 
a(e + e~ Zh°h°) + a{e + e~ -> ZH°) ■ BR(H° h°h ,A°A ) in 2HDM as a function of m H o. 
Right plot: branching ratios of H°. Here m h o = 120 GeV, mjj± = 250 GeV, m^o = 150 GeV, 
?7Ji2 = GeV, tan/3 = 5, and sin a = 0.9. 

H° — > iY ± iJ =F modes, respectively. The cross sections increase by about 20 fb if h°h° and 
A A are produced in e + e~ — > ZH° — > Zh°h° and e + e~ — > — > Zh°h° via resonant 
Higgs-strahlung of H°. 

For the double Higgs-strahlung processes e + e~ — > Zh°h°, e + e~ ZA°A and e + e~ — >• 
Zh°H°, the dominant final states depend on how h°, A and H° decay. In the case where 
h° and A are not heavy (< 125 GeV), h° and A will decay to bb or r + r~, and the final 
states for Zh°h° and ZA°A° would be Z2b2r or Z4r. In the 2HDM, it may be possible 
that, if kinematically allowed, A decays into Zh°. In this case, the final state for ZA°A° 
would be ZA°A° -> 3Zh°h° -> For the e+e" -> Z/i ^ process, the final state will 

be different from the previous one if the Higgs to Higgs decays H° — > h°h° and if — > A A 
are open. As a result, the final state would be Zh°h°h° and Zh°A°A°. 

We would like to stress here that the background study is beyond the scope of this paper. 
However, we point out that in the SM with small cross sections, it has been shown in 
Ref. 10J that for the e + e~ — > ZHH — > Zbbbb process in the SM, the irreducible background 
from electroweak and QCD processes can be suppressed down to manageable levels by using 
kinematics cuts. We expect that in the 2HDM with larger cross sections than in SM, the 
signal can be easily extracted from the background as well. 
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B. Double Higgs-strahlung in the fermiophobic limit 



As commented in Section 2, in 2HDM-I the lightest CP-even Higgs boson h° can be 



fermiophobic 4j, |2l|. This occurs when a = tx/2 (cos a = 0), so that the lightest Higgs 
couplings to all fermions vanish. In this case, the main decay channel of h° (mp o < 2My/) 
is h° — > 77 for m^o < My and h° — > VV* for m^o > my, where V = W, Z {4, 38]. Once the 
WW threshold is crossed, the dominant decay mode becomes h° — > W + W~. 

The LEP Collaboration has already ruled out a fermiophobic Higgs with a mass m^o < 



391 ] . This constraint can be lifted if 



104 GeV and the ZZh° coupling similar to the SM one 
the ZZh° coupling is smaller than the SM one. 

Recently, there is a study devoted to double fermiophobic Higgs bosons production at the 
LHC and ILC 40J. For the ILC, it has been shown that e + e~ — > A°h° followed by the decay 



A —> Zh° can leads to the Zh°h° final state that is similar to our double Higgs-strahlung 
process e + e~ — > Zh°h°. However, the process A°h° — > Zh°h° only depends on the gauge 
couplings while the double Higgs-strahlung process e + e" — > Zh°h° under consideration is 
sensitive to the 2HDM triple Higgs couplings that may enhance the cross section. 

In Fig. O we show the total cross sections of the neutral modes Zh°h° and Zh°H° in 
the fermiophobic limit sin a = 1 for a relatively light CP-even fermiophobic Higgs with 



m^o = 80 GeV. This possibility is not yet ruled out experimentally [39( due to suppressed 
ZZh° coupling. This fermiophobic limit is only relevant for final states having at least one 
fermiophobic Higgs, namely the e + e~ — > Zh°h° and e + e~ — > Zh°H° modes where the h° 
may decays dominantly into two photons. In the left panel of Fig. M we select moderate 
tan/? = 10, and the right panel is for small tan/5 = 1. For both channels, Zh°h° and 
Zh°H°, the cross sections are of the order of few fb and can reach 20 fb (lOfb) for Zh°h° 
in the case of tan/3 = 10 (tan/3 = 1). Due to phase space suppression, the cross section of 
e + e~ — > Zh°H° is smaller than e + e~ — > Zh°h°. 

In the fermiophobic limit of 2HDM-I, one can obtain a very clear signal of 47 + X from 
the Zh°h° mode. For the Zh°H° mode, the signal depends on how H° decays. In 2HDM, 
the heavy Higgs H° can decay to W + W~ or, if kinematically allowed, to h°h°. In the latter 
case, we have ZH°h° — > Zh°h°h° — > Z67 with 67 + X as a distinctive signature. 

To our best knowledge, there is no estimation for the e + e~ — > 47Z backgrounds. We 
expect that such backgrounds should be small. Moreover, requesting four photons in our 
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FIG. 9: The tree level cross sections of the e + e~ — > Zh°hP,Zh H processes as a function of the 
center-of-mass energy in the fermiophobic limit with mi2 = GeV, m^o = 80 GeV, m H o = 160 
GeV, m A o = 300 GeV, m H ± = 250 GeV, and tan/3 = 10 (left) or tan/? = 1 (right). 

signal would be sufficient to kill the backgrounds. 



C. e + e — ► h°h°Z in the decoupling limit 

In the 2HDM, the decoupling limit generally refers to the case when all the scalar masses 
except one formally become infinite and the effective theory is just the SM with one doublet 
(see 15] for recent discussions). In this case, the CP-even h° is the lightest scalar particle 
while the other Higgs particles H°, A and if ± are extremely heavy. Using purely algebraic 
arguments at the tree level, one can derive that the main consequences of the decoupling 
limit are cos(/3 — a) — * 0, and the CP-even h° of the 2HDM and the SM Higgs hsM have 



similar tree-level couplings to gauge bosons and fermions as well 
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431 ] . Obviously, the 



decoupling limit does not rigorously apply to the cases where the particle masses are finite. 



One can consider instead a more realistic scenario, dubbed as the decoupling regime [15], 
where the heavy Higgs particles have masses much larger than the Z boson mass and may 
escape detection in the planned experiments. 

Several studies have been carried out looking for non-decoupling effects in Higgs boson 
decays and Higgs self-interactions. Large loop effects in h° — > 77, h° — * 7Z and h° — > bb 



have been pointed out for the 2HDM 
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421 ] and may provide indirect information on the 



Higgs masses and the involved triple Higgs couplings such as X h o H + H -, \ h o H o H o, A^o^o^o and 
^h°h°h°- The non-decoupling contributions to the triple Higgs self-couplings Xh0h°h° have 
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FIG. 10: a(e + e~ -> Zft°/i°) in units of fb as a function of M$ at ^ = 500 GeV (left) and 
yfs = 800 GeV (right) for different values of m\2- For both plots, m h o = 120 GeV, tan/3 = 2 and 
cos(/3 — a) = 0. 



been investigated in the 2HDM in Ref. 43[, revealing large non-decoupling effects. 

In this section, we will show that the large non-decoupling effects in X h o h o h o modify 
the double Higgs-strahlung e + e~ — > Zh°h° cross section and make it larger than the SM 
expectation. We will focus on the scenario where all the Higgs particles of the 2HDM, except 
for the lightest CP-even Higgs, are heavy and can escape from detection at the first stage 
of next generation colliders. 

It is easy to check that in the decoupling limit (3 — a — > ir/2, the triple Higgs coupling 



^h°h°h° given in Eq. ([5]) reduces to its SM value X^ 0h o = — 3m^ /f- In Ref. |43|], using 



the Feynman diagrammatic method, it has been demonstrated that the one-loop leading 
contributions originated from the heavy Higgs boson loops and the top quark loops to the 
effective A^o/jO^o coupling can be written as 

r2 \ 3 



Ajft 0fc o(2HDM) = ^ 
M 2 N 



H° 



M 2 



in. 



67r 2 m 2 l oV 2 



1 - 



in 



N c mf 
S^m^oV 2 



+ o 



m Jfi/ 

' P 2 m% m% 
m 2 h0 v 2 ' 




(15) 



where M = m^/^sm [3 cos /?, m $ and pi represent the mass of H°, A or H ± and the 
momenta of external Higgs lines, respectively, N c denotes the number of colors, and m t is 
the mass of top quark. We note that in Eq. ( |T5l) m/jO is the renormalized physical mass of 
the lightest CP-even Higgs boson h°. 
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At the tree level, one can show that in the decoupling limit the ZZH° coupling approaches 
and the e + e~ — > Zh°h° amplitude reduces exactly to the SM result. For our calculation of 
the cross section of e + e~ — > Zh°h° in the decoupling limit, we ignore one-loop corrections 
to the ZZh° coupling and replace the h°h°h° coupling by its effective coupling given in 
Eq. (1151) . In fact, it has been shown in Ref. 43j that the non-decoupling effect in ZZh° is 
at a few percent level in the case of cos(/3 — a) = 0. 

In Fig. [TUl we show the cross section of e+e" -> Zh°h° in the SM and 2HDM in its 
decoupling regime as a function of m<j, = m^o = m#o = mjj± . In the calculation of the SM 
cross section, we take into account the one-loop leading contributions originated from top 
quarks given by [43j 



MM^l-JML + O^U, (16) 



For our choice of m^o = 120 GeV, the SM cross section is tiny, « 0.2 fb and 0.16 fb for 
y/s = 500, 800 GeV, respectively. As is clear in Fig. [TQl the 2HDM contributions can enhance 
the cross section by more than one order of magnitude and reach a few fb for small myi and 
large m$. As shown in both panels of Fig. [TU1 we get maximum non- decoupling effect for 
m 12 = 0. It is further amplified with large m$. The plots are cut due to the perturbativity 
and vacuum stability constraints. The perturbativity constraints on \ do not allow m$ 
larger than 850 GeV. In the case of m\i = 500 GeV, the vacuum stability condition does 
not allow m$ to be less than 800 GeV. The sensitivity of cross section to tan j3 is mild, since 
tan/3 only enters the h°h°h° coupling through M 2 = m^ 2 /(sin/3 cos/3). Moreover, tan/3 is 
constrained by perturbativity to be rather moderate. 



IV. SUMMARY 



We have studied the triple Higgs couplings A^ .^ and the double Higgs-strahlung pro- 
duction e + e" — > Z&i§j at linear collider in the framework of general Two Higgs Doublet 
Models. We have quantified the sizes of the triple Higgs couplings compared to the SM 
triple Higgs coupling. We also show that after taking into account the perturbativity and 
vacuum stability constraints on the 2HDM parameters, it is possible to enhance the triple 
Higgs coupling \h,oh°h° U P to 15 times or more than the corresponding SM coupling. If the 
Higgs bosons $ = h°, H°, A and H ± are not too heavy, the double Higgs-strahlung cross 
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sections e + e~ — > Z$i&j can be substantial, at the level of a few hundred fb, and provide 
some information on the triple Higgs couplings. 

We have also studied the double Higgs-strahlung processes e + e" — > Zh°h° and e + e~ — > 
ZH°h° in the fermiophobic limit of 2HDM-I where the collider signature can be very dis- 
tinctive with the final state Z2^y, ZAry, or Z67 if H° — > h°h°. We also analyze the double 
Higgs-strahlung process e + e~ — * Zh°h° in the decoupling limit where h° mimics the SM 
Higgs boson. It is shown that in this limit, the cross section can be enhanced by about two 
orders of magnitude, which is much larger than the MSSM enhancement. Observations of 
such large cross sections would definitely indicate that the Higgs sector is 2HDM-like. 
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